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INTRODUCTION 

Traumatic  brain  injury  (TBI)  is  a  major  casualty  identified  among  veterans  deployed  to  the 
Persian  Gulf  region  in  support  of  Operation  Enduring  Freedom  (OEF)  or  Operation  Iraqi 
Freedom  (OIF).  TBI  is  caused  by  one  or  more  concussive  insults  to  the  head  or  a  penetrating 
head  injury  that  disrupts  the  normal  function  of  the  brain,  leading  to  either  transient  or  chronic 
impairments  in  physical,  cognitive,  emotional  and  behavioral  functions.  TBI  in  OEF/OIF 
veterans  is  largely  the  result  of  concussive  injuries  from  blast-producing  weaponry.  Overall, 
veterans  have  lower  attention  scores,  although  it  is  not  yet  known  if  mild  TBI  might  contribute  to 
this  observation.  Nevertheless,  veterans  from  prior  conflicts  exposed  to  blasts  have  shown 
evidence  of  mild  TBI  and  attention  difficulties  when  compared  to  similar  veterans  without  blast 
exposure.  Early  diagnosis  of  chronic  TBI  is  important  for  preventing  further  progression  of 
symptoms  that  can  disrupt  a  veteran’s  life  upon  return  from  service  overseas.  Mild  TBI  can  be 
difficult  to  diagnose  and,  when  coupled  with  psychological  illness,  can  be  either  misdiagnosed  or 
missed  altogether.  Traditionally,  physicians  and  scientists  have  viewed  and  interpreted  diseases 
at  the  ‘visuaF  clinical  level.  With  the  advent  of  genomics  and  proteomics  technologies, 
personalized  medicine  offers  the  promise  and  potential  of  uncovering  the  largely  ‘unseen’  details 
of  disease  causality,  onset,  and  progression.  As  we  have  proposed,  these  studies  were  designed  to 
be  conducted  in  collaboration  with  The  War-related  Illness  and  Injury  Study  Center  (WRIISC), 
Department  of  Veteran  Affairs,  New  Jersey  Health  Care  System  (DVANJHCS),  East  Orange, 

NJ,  and  were  designed  to  identify  genomic-microRNA  fingerprints  from  clinically  accessible 
blood  cell  components  as  independent  biological  indexes  that  will  allow  us  to  identify  unique 
molecular  indices  of  Persistent  Postconcussive  Syndrome  as  well  as  well  as  other  significant 
injury-related  factors  associated  with  mild  TBI  in  OEF/OIF  veterans. 

BODY 

The  study  was  designed  to  identify  clinically  accessible  molecular  biomarkers  of  TBI  prior  to 
definitive  clinical  diagnosis  using  high  throughput  microRNA  technology.  Specifically,  these 
studies  were  meant  to  identify,  characterize,  and  validate  microRNA  biomarker  species  whose 
content  in  peripheral  blood  mononuclear  cells  (PBMC)  could  help  to  distinguish  TBI  injury 
cases  within  a  veteran  population  following  deployment  in  support  of  OEF  or  OIF.  The  overall 
study  is  separated  into  a  Biomarker  Discovery  study  to  identify  candidate  microRNA 
biomarkers,  and  a  Biomarker  Validation  study  to  validate  the  sensitivity  and  specificity  of 
microRNA  biomarkers,  either  individually  or  as  panels  of  multiple  biomarkers,  to  correctly 
identify  TBI  and  control  cases  in  an  independent  cohort  of  TBI  and  non-TBI  veterans. 

The  WRIISC  at  East  Orange,  New  Jersey  received  final  Department  of  Defense  (DoD)  IRB 
approval  to  commence  recruitment  on  July  29,  2009.  As  was  stated  in  our  Year  3  Annual 
Progress  Report,  which  was  submitted  in  May  20 1 1 ,  we  initiated  volunteer  recruitment  on 
October  3 1 ,  2009.  Our  inclusion  criteria  for  recruitment  were  male  and  females  1 8-75  years  with 
or  without  a  history  of  TBI,  who  have  completed  a  clinical  evaluation  at  the  East  Orange,  New 
Jersey  WRIISC.  Participants  were  classified  as  having  a  history  of  mTBI  if  they  positively 
endorsed  at  least  one  of  4  items  on  the  Veteran  traumatic  brain  injury  screening  tool  (VAT- 
BIST)  (Donnelly  et  ah,  2011)  and  had  a  score  at  least  one  standard  deviation  below  the  norm  for 
age  and  education  on  the  Repeatable  Battery  for  Neuropsychological  Testing  (RBANS) 
(Randolph,  1998).  Classification  criteria  for  Control  cases  are  DVBIC  confirmation  of  no  injury 
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to  the  head  and  an  RBANS  score  less  than  one  standard  deviation  below  the  nonn.  Cases  with 
inter-current  infections  of  inflammatory-related  conditions  were  excluded. 

Starting  on  October  31,  2009,  when  recruitment  was  initiated,  and  continuing  through  the  time 
when  we  submitted  our  Year  3  Annual  Report  in  May,  201 1,  we  recruited  and  collected  PBMC 
specimens  from  97  cases,  including  22  TBI  and  75  non  TBI  cases  (Table  I).  The  proportion  of 
veterans  classified  as  mild  TBI  in  our  recruited  cohort  is  consistent  with,  and  even  slightly 
above,  previous  prevalence  estimates  of  12%  (Schneiderman  et  al,  2008)  reported  in  a  cross- 
sectional  survey  of  2,235  active  duty,  guard  and  reserve  OEF/OIF  veterans. 


Table  I.  The  number  of  mTBI 
and  non-TBI  control  veterans 
recruited  during  Years  1,  2  and 

3.  A  total  of  97  OEF/OIF  veteran 
cases  were  recruited  in  Years  1-3, 
among  which  22  (22.6%)  were 
mTBI  cases.  Also  shown  are  the 
numbers  of  cases  used  for 
Biomarker  Discovery  studies  as 
well  as  the  number  of  cases  available  for  Biomarker  Validation  studies 


mTBI 

Non-TBI 

control 

Total 

Recruitment:  Years  1-3 

22 

75 

97 

Cases  used  for 
Biomarker  Discovery  studies 

11 

16 

27 

Cases  available  for  Biomarker 

Validation  studies  at  the  end  of  Year  3 

11 

59 

70 

From  these  97  cases  recruited  during  Years  1-3,  we  used  1 1  TBI  and  16  non-TBI  cases  for  an 
Interim  Biomarker  Discovery  study.  As  RNA  preparations  from  some  of  the  cases  did  not  meet 
quality  control  criteria,  our  final  Biomarker  Discovery  studies  were  conducted  during  Year  3 
using  9  TBI  and  9  non  TBI  cases.  At  the  time  when  we  submitted  our  Year  3  Annual  Report,  we 
had  1 1  TBI  and  59  non-TBI  (a  total  of  70)  cases  that  were  not  used  in  our  Biomarker  Discovery 
studies  and  therefore  were  available  for  Biomarker  Validation  studies  (Table  I). 

Details  of  our  interim  microRNA  Biomarker  Discovery  study  were  provided  in  our  Year  3 
Annual  Progress  Report  that  we  submitted  in  May,  2011.  Key  outcomes  from  our  Biomarker 
Discovery  studies  are  as  follows:  Using  a  high  throughput  array  chip  platform,  we  initially 
identified  18  candidate  small  RNA  TBI  biomarkers  from  our  Veteran  study  cohort.  Using 
independent  Q-PCR  assays,  we  confirmed  that  13  of  the  18  candidate  small  RNA  biomarkers 
are,  indeed,  differentially  regulated  in  the  PBMC  of  TBI  compared  to  non-TBI  veteran  cases. 
The  13  confirmed  small  RNA  biomarkers  include  12  small  nucleolar  RNA  (ACA48, 
ENSG19941 1,  HBII-239,  HBII-289,  U15B,  U27,  U35A,  U55,  U56,  U58B,  U83A,  U91)  and  1 
miRNA  (Has-miR-671-5p).  Each  of  the  13  confirmed  small  RNA  biomarkers  are  found  in 
significantly  lower  levels  in  PBMC  specimens  from  TBI  compared  to  non-TBI  control  veteran 
cases.  Using  unsupervised  clustering,  we  found  that  the  differential  regulation  of  these  small 
RNA  biomarkers  in  clinically  accessible  blood  cells  is  capable  of  distinguishing  TBI  from  non- 
TBI  OEF/OIF  Veteran  cases.  Moreover,  we  identified  a  3-biomarker  panel  capable  of 
distinguishing  TBI  from  non-TBI  veteran  cases  with  89%  accuracy,  82%  selectivity  and  78% 
specificity. 
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Based  on  outcomes  from  our  Biomarker  Discovery  studies  identifying  and  confirming  13  small 
RNA  biomarkers  from  our  Biomarker  Discovery  cohort,  we  are  continuing  with  our  Biomarker 
Validation  studies  to  test  the  value  of  these  biomarkers  in  tenns  of  correctly  diagnosing  TBI  and 
non-TBI  cases  in  a  new  independent  cohort  of  Veterans  as  we  proposed.  Based  on  power 
analysis  using  the  real-time  quantitative  PCR  (Q-PCR)  data  from  the  Biomarker  Discovery  study 
cohort,  we  calculated  that  a  group  size  of  23  TBI  cases  and  an  equal  number  of  non-TBI  cases 
will  provide  a  power  of  0.99  (alpha  level  =  0.05)  for  detecting  differences  in  each  of  the  13  small 
RNA  biomarkers  in  TBI  compared  to  non-TBI  OEF/OIF  cases.  Since  the  end  of  Year  3  we  had 
only  1 1  TBI  and  non-TBI  cases  that  were  available  for  Biomarker  Discovery  studies  (Tables  I 
and  II),  we  applied  for  a  one  year  no-cost  extension  so  that  we  could  continue  to  recruit  OIF/OEF 
veteran  cases  and  complete  our  proposed  Biomarker  Validation  studies.  We  received  approval 
for  a  second  no-cost  extension  on  March  5th,  2012,  allowing  us  to  continue  with  the  proposed 
studies  to  April  30,  2013. 

Our  collaborators  at  The  War  Related  Illness  and  Injury  Study  Center  in  the  East  Orange  New 
Jersey  Health  Care  System  (WRIISC/NJHCS)  continued  to  recruit  OEF/OIF  veterans.  Following 
the  submission  of  our  Year  3  Annual  Report  in  May  201 1,  we  recruited  an  additional  13  new 
cases  (6  TBI  and  7  non-TBI  cases)  for  a  total  of  83  cases,  comprised  of  17  TBI  and  66  non-TBI 
control  cases,  which  are  available  for  Biomarker  Validation  studies  (Table  II).  No  problems  were 
encountered  with  volunteer  recruitment  or  blood  collection  and  processing. 


Table  II.  Cases  currently 
available  for  Biomarker 
Validation  studies.  Shown 
are  the  numbers  ofmTBI 
and  control  OEF/OIF 
cases  recruited  in  Years  1- 
3  and  in  Year  4,  and  the 
number  of  TBI  and  control 
cases  that  are  currently 
available  for  TBI 

Biomarker  Validation  studies. 

As  we  have  noted  above,  power  analysis  calculated  that  Biomarker  Validation  studies  conducted 
using  PBMC  specimens  from  23  TBI  and  23  age-  and  gender-matched  non-TBI  cases  will 
provide  a  power  of  0.99  (alpha  level  =  0.05)  to  detect  differences  in  the  levels  of  the  small  RNA 
biomarkers  in  TBI  compared  to  non-TBI  cases.  Nonetheless,  based  on  power  analysis,  our 
presently  available  17  TBI  cases  and  66  control  cases  will  provide  a  power  of  0.95  (alpha  level  = 
0.05)  for  detecting  differences  in  the  contents  of  small  RNA  TBI  biomarkers  in  PBMC 
specimens  from  TBI  compared  to  non-TBI  cases.  The  WRIISC/NJHCS  has  sent  samples  from  all 
currently  available  83  veteran  cases  by  commercial  courier  to  Dr.  Pasinetti’s  group  at  the  JJ 
Peters  VA  Medical  Center  for  Biomarker  Validation  studies.  Q-PCR  analysis  of  the  83  cases 
Biomarker  Validation  study  cohort  is  ongoing.  We  also  are  continuing  to  recruit  additional  cases 
for  our  Biomarker  Validation  studies.  As  we  have  done  in  the  past  year,  we  anticipate  we  will 
recruit  additional  13  new  cases  in  the  upcoming  year.  Based  on  our  estimation  that  25%  of  the 


mTBI 

Non-TBI 

control 

Total 

Cases  available  for  Biomarker 

Validation  studies  at  the  end  of  Year  3 

11 

59 

70 

New  cases  recruited  in  Year  4 

6 

7 

13 

Cases  currently  available  for 

Biomarker  Validation  studies 

17 

66 

83 
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recruited  cases  will  be  identified  as  TBI,  we  this  will  increase  the  number  of  TBI  cases  to  20  for 
Biomarker  Validation  studies  (Table  III),  which  will  increased  our  likelihood  (to  97%)  of 
detecting  differences  in  the  regulation  of  small  TBI  biomarkers  in  TBI  vs.  non- TBI  cases. 


mTBI 

Non-TBI 

control 

Total 

Cases  currently  available  for 
Biomarker  Validation  studies 

17 

66 

83 

Anticipated  13  new  cases 
to  be  recruited  in  the  upcoming  year 

3 

10 

13 

Total  number  of  cases  that  will 

be  used  for  Biomarker  Validation  studies 

20 

76 

96 

Table  III:  Addition  of  13  new 
recruits  will  provide  a  total  of 
20  TBI  cases  for  Biomarker 
Validation  studies. 


In  Year  4,  we  also  prepared  a  manuscript  entitled  “ Select  non-coding  RNA  in  blood  components 
provide  novel  clinically  accessible  biological  surrogates  for  improved  identification  of  traumatic 
brain  injury  in  OEF/OIF  Veterans ”  based  on  observations  from  our  TBI  Biomarker  Discovery 
studies.  Our  manuscript  was  recently  accepted  for  publication  by  the  American  Journal  of 
Neurodegenerative  Disease.  In  our  manuscript,  we  noted  that  the  majority  of  TBI  cases  in  our 
study  cohort  are  co-morbid  with  PTSD  based  on  a  PTSD  diagnosis  criterion  of  having  a  score  of 
50  or  more  in  the  PTSD  Checklist  -  Civilian  Version  (Harris  et  ah,  2008),  and  that  our  non-TBI 
control  cases  were  selected  to  match  for  PTSD  diagnosis.  Thus,  our  identified  panel  of  13  small 
RNA  biomarkers  likely  represents  biological  indices  selective  for  TBI.  Moreover,  we  noted  that 
TBI  cases  in  our  biomarker  discovery  studies  were  recruited  after  an  average  interval  of  3.9  years 
following  their  last  deployment  (deployment-to-recruitment  interval  ranging  from  0.7  to  10.2 
years,  with  a  median  interval  of  3.4  years).  Thus,  changes  in  the  regulation  of  these  small  RNA 
TBI  biomarkers  that  we  observed  are  not  acute  TBI  responses,  but  likely  represent  long-term 
pathophysiological  consequences  subsequent  to  TBI.  Lastly,  we  observed  that  genes  relevant  to 
neural  circuits,  synapses  and  neural  plasticity  processes  are  also  expressed  in  circulating  blood 
cells,  such  as  PBMC.  Thus,  significant  down-regulation  of  select  small  RNA  biomarkers  we 
observed  in  PBMC  specimens  from  our  Veteran  TBI  population  long  after  their  deployment 
might  reflect  long-temi  molecular  alterations  in  the  central  nervous  system  contributing  to  the 
onset  and  progression  of  clinical  TBI  phenotypes.  We  attached  a  copy  of  our  TBI  biomarker 
manuscript  as  an  appendix  (Appendix  I)  accompanying  our  Year  4  Annual  Progress  Report. 


KEY  RESEARCH  ACCOMPLISHMENTS 

1 .  We  conducted  Biomarker  Discovery  studies  using  a  high-throughput  Array  chip  platfonn 
and  identified  18  candidate  TBI  small  RNA  biomarkers. 

2.  Using  independent  Q-PCR  assays,  we  confirmed  that  13  of  these  candidate  small  RNA 
biomarker  species  are,  indeed,  significantly  down-regulated  in  PBMC  of  TBI  compared 
to  non-TBI  control  veteran  cases. 

3.  Using  unsupervised  clustering,  we  found  that  the  differential  regulation  of  these  small 
RNA  biomarkers  in  clinically  accessible  blood  cells  is  capable  of  differentiating  between 
TBI  and  non-TBI  OEF/OIF  veteran  cases. 
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4.  We  identified  a  3-biomarker  panel  capable  of  distinguishing  TBI  from  non-TBI  control 
veteran  cases  with  89%  accuracy  82%  selectivity  and  78%  specificity. 

5.  The  majority  of  TBI  cases  in  our  biomarker  study  were  co-morbid  with  PTSD  and  our 
non-TBI  control  cases  were  selected  to  match  for  PTSD  diagnosis.  Thus,  our  identified 
panel  of  13  small  RNA  biomarkers  likely  represents  biological  indices  selective  for  TBI. 

6.  We  have  recruited  OEF/OIF  Veteran  volunteers.  We  now  have  a  total  of  17  TBI  and  66 
non-TBI  control  cases  for  our  proposed  Biomarker  Validation  studies.  We  calculated  that 
the  number  of  currently  available  cases  will  provide  a  power  of  0.95  (alpha  level  =  0.05) 
for  detecting  differences  in  the  contents  of  small  RNA  TBI  biomarkers  in  PBMC 
specimens  from  TBI  compared  to  non-TBI  cases. 

7.  We  have  published  one  manuscript  on  the  identification  and  characterization  of  small 
RNA  TBI  biomarkers  from  clinically  accessible  PBMC. 


REPORTABLE  OUTCOMES 

N/A 


CONCLUSION 

Our  Biomarker  Discovery  studies  to  date  have  identified  a  panel  of  small  RNA  TBI  biomarkers 
whose  contents  in  clinically  accessible  blood  cells  provides  a  criterion  for  the  correct  diagnosis 
of  TBI  from  non-TBI  veteran  cases  with  high  accuracy,  specificity,  selectivity.  The  majority  of 
TBI  cases  in  our  biomarker  study  were  co-morbid  with  PTSD  and  our  non-TBI  control  cases 
were  selected  to  match  for  PTSD  diagnosis.  Thus,  our  identified  panel  of  13  small  RNA 
biomarkers  likely  represents  biological  indices  selective  for  TBI.  In  Year  4,  we  recruited  an 
additional  13  new  cases  (6  TBI  and  7  non-TBI  cases)  for  a  total  of  83  cases,  comprised  of  17  TBI 
and  66  non-TBI  control  cases,  which  are  available  as  an  independent  study  cohort  for  Biomarker 
Validation  studies.  We  calculated  that  the  currently  available  83  cases  are  sufficiently  powered 
for  our  proposed  Biomarker  Validation  studies.  In  ongoing  studies  we  are  analyzing  the  PBMC 
specimens  from  the  Biomarker  Validation  study  cohort  using  Q-PCR  to  validate  individual  or 
combination  of  small  RNA  TBI  biomarkers  as  surrogate  biological  indices  of  TBI. 

Outcomes  from  our  combined  Biomarker  Discovery  and  Validation  studies  will  lead  to  improved 
TBI  detection  in  OEF/OIF  veterans,  more  sensitive  outcome  measurements  for  future  clinical 
trials,  a  better  understanding  of  the  biological  mechanisms  underlying  concussive  TBI,  and 
insights  into  novel  therapeutic  targets  for  OEF/OIF  veterans  with  TBI. 
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Abstract:  This  study  was  designed  to  identify  clinically  accessible  molecular  biomarkers  of  mild  traumatic  brain  injury 
(mTBI)  that  could  be  used  to  help  identify  returning  Operation  Iraqi  Freedom  (OIF)  and  Operation  Enduring  Freedom 
(OEF)  Veterans  who  are  suffering  from  the  effects  of  mTBI.  While  analyzing  the  expression  profile  of  small  non-coding 
RNAs  in  peripheral  blood  mononuclear  cells  (PBMCs)  from  an  OEF/OIF  veteran  study  cohort  using  a  high  throughput 
array  chip  platform,  we  identified  18  candidate  small  non-coding  RNA  biomarkers  that  are  differentially  regulated  in 
PBMCs  of  mTBI  compared  to  non-TBI  control  cases.  Independent  quantitative  real-time  polymerase  chain  reaction 
assays  confirmed  that  13  of  these  candidate  small  RNA  biomarker  species  are,  indeed,  significantly  down-regulated 
in  PBMCs  of  mTBI  compared  to  non-TBI  control  veteran  cases.  Based  on  unsupervised  clustering  analysis,  we  identi¬ 
fied  a  3-biomarker  panel  which  was  most  able  to  distinguish  mTBI  from  non-TBI  control  veteran  cases  with  high  accu¬ 
racy,  selectivity  and  specificity.  The  majority  of  mTBI  cases  in  our  biomarker  study  were  co-morbid  with  Post- 
Traumatic  Stress  Disorder  (PTSD),  and  thus  our  non-TBI  control  cases  were  selected  to  match  PTSD  diagnoses. 
Therefore,  our  identified  panel  of  3  small  RNA  biomarkers  likely  represents  a  biological  index  selective  for  mTBI.  Out¬ 
comes  from  our  studies  suggest  that  additional  applications  of  the  clinically  accessible  small  non-coding  RNA  bio¬ 
markers  to  current  diagnostic  criteria  may  lead  to  improved  mTBI  detection  and  more  sensitive  outcome  measures 
for  clinical  trials.  Future  studies  exploring  the  physiological  relevance  of  mTBI  biomarkers  will  also  provide  a  better 
understanding  of  the  biological  mechanisms  underlying  mTBI  and  insights  into  novel  therapeutic  targets  for  mTBI. 

Keywords:  Mild  traumatic  brain  injury  (mTBI),  biomarkers,  microRNA  (miRNA),  post-traumatic  stress  disorder  (PTSD) 


Introduction 

Traumatic  brain  injury  (TBI)  is  a  condition  often 
identified  among  Veterans  deployed  to  the  Per¬ 
sian  Gulf  region  in  support  of  Operation  Endur¬ 
ing  Freedom  (OEF)  or  Operation  Iraqi  Freedom 
(OIF).  TBI  is  caused  by  one  or  more  concussive 
insults  to  the  head  or  a  penetrating  head  injury 
that  disrupts  the  normal  functions  of  the  brain, 
leading  to  either  transient  or  chronic  impair¬ 
ments  in  physical,  cognitive,  emotional  and  be¬ 
havioral  functions  [1-6].  In  OEF/OIF  Veterans, 
TBI  is  largely  the  result  of  concussive  injuries 
from  blast-producing  weaponry  [7].  Veterans 
exposed  to  blasts  from  prior  conflicts  have 


shown  evidence  of  mild  TBI  (mTBI)  and  atten¬ 
tion  difficulties  when  compared  to  similar  Veter¬ 
ans  without  blast  exposure.  Mild  TBI  can  be  dif¬ 
ficult  to  diagnose  and,  when  coupled  with  psy¬ 
chological  illness,  can  be  either  misdiagnosed 
or  missed  altogether.  Traditionally,  physicians 
and  scientists  have  viewed  and  interpreted  dis¬ 
eases  at  the  ‘visual’  clinical  level.  However,  with 
the  advent  of  genomics  and  proteomics  tech¬ 
nologies,  personalized  medicine  offers  the 
promise  and  potential  of  uncovering  the  largely 
‘unseen’  details  of  disease  causality,  onset,  and 
progression. 

New  evidence  has  highlighted  defects  in  neural 
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circuits  and  synapses,  and  the  plastic  processes 
controlling  these  functions,  in  TBI  [8-13].  While 
gene  products  relevant  to  these  processes  are 
expressed  in  the  brain,  some  of  these  genes  are 
also  expressed  in  circulating  blood  cells,  such 
as  peripheral  blood  mononuclear  cells  (PBMCs) 
[14-17].  Consistent  with  this,  recent  studies 
have  illustrated  that  PBMC-associated 
biomarkers  may  provide  insights  into  the 
pathogenesis  of  neurological  disorders  such  as 
Alzheimer’s  disease  and  may  be  used  to 
monitor  disease  diagnosis  and  progression 
[18,19].  Thus,  PBMCs  may  also  provide  an  ideal 
clinically  accessible  “window”  into  the  brain, 
reflecting  molecular  alterations  following  TBI 
which  might  contribute  to  the  onset  and 
progression  of  clinical  TBI  phenotypes. 

Small  non-coding  RNAs,  including  microRNA 
(miRNA)  and  small  nucleolar  RNA  (snoRNA),  are 
increasingly  recognized  for  their  roles  in  the 
regulation  of  cellular  processes  in  health  and 
disease  [20].  Select  small  non-coding  RNAs, 
particularly  miRNA,  have  been  implicated  in 
neurological  disorders  [21-23].  It  is  possible 
that  miRNA  and  other  small  non-coding  RNAs 
might  contribute  to  the  onset  and/or  progres¬ 
sion  of  clinical  complications  following  TBI  [24]. 
Exploring  the  feasibility  of  identifying  clinically 
accessible  TBI  biomarkers,  we  identified  select 
small  non-coding  RNA  fingerprints  from  clinically 
accessible  PBMCs  that  may  be  used  as  inde¬ 
pendent  biological  indexes  of  mTBI  in  OEF/OIF 
Veterans.  Outcomes  from  our  studies  suggest 
that  additional  applications  of  the  clinically  ac¬ 
cessible  small  non-coding  RNA  biomarkers  to 
current  diagnostic  criteria  may  lead  to  improved 
mTBI  detection  and  more  sensitive  outcome 
measures  for  clinical  trials. 

Materials  and  methods 

18  OIF  and  OEF  Veterans  (9  mTBI  and  9  non- 
mTBI  control  cases)  were  recruited  by  The  War 
Related  Illness  and  Injury  Study  Center 
(WRIISC),  Department  of  Veterans  Affairs,  New 
Jersey  Health  Care  System  (DVANJHCS),  East 
Orange,  NJ.  Male  and  female  participants  were 
included  if  they  were  between  18-75  years  of 
age  and  completed  a  clinical  evaluation  at  the 
New  Jersey  WRIISC.  Participants  were  included 
regardless  of  their  mTBI  history.  Cases  with  inter 
-current  infections  or  inflammatory-related  con¬ 
ditions  were  excluded.  Participants  were  classi¬ 
fied  as  having  a  history  of  mTBI  if  they  posi- 
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tively  endorsed  at  least  one  of  4  items  on  the 
Veteran  traumatic  brain  injury  screening  tool 
(VAT-BIST)  [25]  and  had  a  score  at  least  one 
standard  deviation  below  the  norm  for  age  and 
education  on  the  Repeatable  Battery  for  Neuro¬ 
psychological  Testing  (RBANS)  [26]  .  Classifica¬ 
tion  criteria  for  Control  cases  included  a  nega¬ 
tive  VAT-BIST  score  and  a  RBANS  score  less 
than  one  standard  deviation  below  the  norm. 

Demographic  information  for  individual  mTBI 
and  non-mTBI  cases  is  presented  in  Table  1. 
The  average  age  of  mTBI  and  non-mTBI  cases 
used  in  our  interim  Biomarker  Discovery  study 
was,  respectively,  31.6±7.0  and  29.8+8.2 
years.  The  interval  between  Veterans’  last  de¬ 
ployment  and  recruitment  into  this  study  was 
3.9±2.7  and  2.6±2.1  years  for  the  mTBI  and 
non-mTBI  group,  respectively.  The  mTBI  group 
had  an  average  of  13.3H.3  years  of  education 
and  the  non-mTBI  group  had  an  average 
13.0±2.4  years  of  education.  There  was  no  sig¬ 
nificant  difference  in  age,  deployment  interval 
or  years  of  education  between  the  mTBI  and  the 
non-mTBI  groups  (t-test  assessments  of  mTBI 
versus  non-mTBI  groups:  p-values  of  0.59  for 
age,  0.30  for  deployment  interval,  and  0.72  for 
duration  of  education).  The  proportion  of  males 
in  the  mTBI  and  the  non-mTBI  group  was,  re¬ 
spectively,  78%  and  67%.  Lastly,  89%  of  the 
mTBI  veteran  cases  used  in  our  study  were  co- 
morbid  with  post-traumatic  stress  disorder 
(PTSD),  based  on  a  PTSD  diagnosis  criterion  of 
having  a  score  of  50  or  more  in  the  PTSD 
Checklist  -  Civilian  Version  [27].  Thus  non-mTBI 
cases  were  selected  to  match  for  PTSD,  with 
78%  of  cases  in  the  non-mTBI  control  group 
diagnosed  with  PTSD. 

PBMC  isolation 

Blood  specimens  were  collected  by  venipunc¬ 
ture  and  drawn  into  BD  Vacutainer  CPT  Cell 
Preparation  Tubes.  PBMCs  were  isolated  from 
freshly  collected  blood  specimens  following 
manufacturer's  instructions  (Becton,  Dickinson 
and  Company)  and  were  stored  at  -80  °C  until 
use. 

RNA  preparation  and  high  throughput  analysis 
of  small  non-coding  RNAs 

Total  RNA  was  isolated  from  approximately  10- 
50  mg  of  PBMCs  using  RNA  STAT-60  according 
to  the  manufacturer's  instructions  (Tel-Test, 
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Table  1.  Demographic  characteristics  of  mTBI  and  non-mTBI  control  cases  we  used  in  our  interim  Biomarker  Discov¬ 
ery  study.  Mild  TBI  classification  is  based  on  positive  endorsement  of  the  VA  traumatic  brain  injury  screen  (VAT  BIS) 
and  a  score  of  at  least  one  standard  deviation  below  the  norm  for  age  and  education  on  the  Repeatable  Battery  for 
Neuropsychological  Testing  (RBANS).  Non-mTBI  case  classification  is  based  on  negative  endorsement  of  VAT-BIS 
and  a  RBANS  score  of  less  than  one  standard  deviation  below  the  norm.  PTSD  diagnosis  is  based  on  a  score  of  50 
or  more  on  the  PTSD  Checklist  -  Civilian  Version.  Average  age:  mTBI  group,  31.6±7.0  yrs;  non-mTBI  group, 
29.8±8.2  yrs.  Interval  between  their  last  deployment  and  recruitment  into  this  study:  mTBI  group,  3.9±2.7  yrs;  non- 
mTBI  group  2.6±2.1  yrs.  Average  duration  of  education:  mTBI  group,  13.3±1.3;  non-mTBI  group,  13.0±2.4  yrs.  Per¬ 
centage  of  males:  mTBI  group,  78%;  non-mTBI  group,  67%.  Percent  of  cases  with  co-morbid  PTSD:  mTBI  group,  89%; 
non-mTBI  group,  78%. 


Case 

mTBI/Ctl 

Age 

Gender 

Ethnicity 

Interval  (yrs)  since 
last  deployment 

Education 

(yrs) 

Comorbidity 

PTSD 

31529 

mTBI 

38 

Male 

Black,  non- 
Hispanic 

3.0 

14 

Yes 

33297 

mTBI 

41 

Male 

Native  Ameri-  4.3 

can 

12 

Yes 

33825 

mTBI 

31 

Male 

Black,  non- 
Hispanic 

3.4 

16 

Yes 

33828 

mTBI 

42 

Male 

Black,  non- 
Hispanic 

0.7 

14 

Yes 

33888 

mTBI 

27 

Female 

White  His¬ 
panic 

4.5 

14 

Yes 

33931 

mTBI 

23 

Male 

White  His¬ 
panic 

1.2 

12 

Yes 

33947 

mTBI 

25 

Female 

Black,  non- 
Hispanic 

2.8 

13 

No 

33881 

mTBI 

32 

Male 

White  His¬ 
panic 

10.2 

13 

Yes 

33811 

mTBI 

27 

Male 

Black,  non- 
Hispanic 

4.8 

12 

Yes 

31705 

Non  TBI  Ctl 

27 

Male 

Black,  non- 
Hispanic 

3.4 

12 

No 

33565 

Non  TBI  Ctl 

25 

Male 

White  His¬ 
panic 

1.1 

12 

Yes 

33578 

Non  TBI  Ctl 

30 

Female 

White  His¬ 
panic 

4.3 

16 

Yes 

33596 

Non  TBI  Ctl 

35 

Male 

White  His¬ 
panic 

3.8 

16 

Yes 

33598 

Non  TBI  Ctl 

49 

Male 

White  His¬ 
panic 

0.7 

9 

Yes 

33821 

Non  TBI  Ctl 

26 

Female 

White  His¬ 
panic 

0.2 

16 

No 

33834 

Non  TBI  Ctl 

24 

Male 

White  His¬ 
panic 

2.8 

12 

Yes 

33913 

Non  TBI  Ctl 

22 

Female 

Black,  non- 
Hispanic 

6.5 

12 

Yes 

33930 

Non  TBI  Ctl 

30 

Male 

Black,  non- 
Hispanic 

1.0 

12 

Yes 

Friendswood, 

TX,  USA). 

Immediately  prior 

to 

Total  RNA  was  directly  modified 

and  labeled 

RNA  labeling, 

the  purity  and 

concentration 

of 

using  the  FlashTag™ 

HSR  Biotin 

RNA  Labeling 

RNA  samples  were  determined  from  OD260/280 
readings  using  a  dual  beam  UV  spectropho¬ 
tometer  and  RNA  integrity  was  determined  by 
capillary  electrophoresis  using  the  RNA  6000 
Nano  Lab-on-a-Chip  kit  and  the  Bioanalyzer 
2100  (Agilent  Technologies,  Santa  Clara,  CA, 
USA)  as  per  the  manufacturer’s  instructions. 


Kit  according  to  the  manufacturer’s  instructions 
(Genisphere,  Hatfield,  PA).  Verification  of  biotin 
labeling  was  obtained  by  an  enzyme-linked  oli- 
goabsorbant  assay  (ELOSA)  using  Immobilizer™ 
Amino  -  8  well  strips  (Nunc/Thermo  Fisher  Sci¬ 
entific,  Rochester,  NY,  USA)  according  to  instruc¬ 
tions  supplied  by  Genisphere.  Labeled  cRNA 
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(1.0  (jg)  was  hybridized  for  16hr  at  48°C  to  Affy- 
metrix  microRNA  vl.O  arrays  (Affymetrix,  Santa 
Clara,  CA,  USA),  which  contain  probe  sets  for 
1,500  small  non-coding  RNAs,  including  mi¬ 
croRNA  (miRNA),  small  nucleolar  RNA  (snoRNA), 
small  Cajal  body-specific  RNA  (scaRNA),  and 
5.8S  ribosomal  RNA  (rRNA).  Array  content  was 
derived  from  the  Sanger  miRBase  miRNA  data¬ 
base  vll  (April  15,  2008,  http:// 

microrna.sanger.ac.uk),  snoRNABase 
(www.snorna.biotoul.fr)  and  the  Ensembl  Ar¬ 
chive  (www.ensembl.org).  Arrays  were  washed 
and  stained  on  a  Fluidics  Station  450 
(Affymetrix)  according  to  the  manufacturer’s 
recommended  procedures.  The  arrays  were 
stained  with  phycoerythrein-conjugated  strepta- 
vidin  (Invitrogen/Life  Technologies,  Carlsbad, 
CA,  USA)  and  the  fluorescence  intensities  deter¬ 
mined  using  a  GCS  3000  7G  high-resolution 
confocal  laser  scanner  and  AGCC  software 
(Affymetrix).  The  scanned  images  were  analyzed 
with  the  miRNA  QC  tool  (Affymetrix)  using  RMA 
global  background  correction,  quantile  normali¬ 
zation  and  median  polish  summarization  to  gen¬ 
erate  quantified  data  (as  recommended  by 
Genisphere).  Quality  control  metrics  for  arrays 
included  normalized  signal  values  >  1000  for 
five  spike-in  control  oligo  probe  sets 
(Genisphere). 

Small  RNA  probe  sets  exhibiting  significant  dif¬ 
ferential  expression  (SDE)  were  identified  using 
the  following  steps  in  GeneMaths  XT  (Applied 
Maths,  Austin  TX):  1)  Probed  sets  with  array 
detection  p-values  <  0.05  for  all  samples  in  at 
least  one  experimental  group  were  selected  for 
further  analysis,  2)  Performed  Discriminant 
Analysis  (DA)  and  determined  the  largest  per¬ 
centage  of  remaining  probe  sets  that  permitted 
correct  group  assignment  of  samples  in  unsu¬ 
pervised  hierarchical  clustering  by  the  Un¬ 
weighted  Pair-Group  Method  using  Arithmetic 
averages  (UPGMA)  based  on  cosine  correlation 
of  row  mean  centered  log2  signal  values;  this 
was  the  top  50%-tile,  3)  In  the  DA  top  50%-tile, 
selected  probe  sets  with  absolute  signal  log2 
fold  changes  >  1.0  and  independent  t-test  p- 
values  <  0.05  adjusted  for  multiple  testing  error 
by  the  Benjamini-Hochberg  false-discovery  rate 
(FDR)  correction  method  [28].  Unsupervised 
hierarchical  clustering  of  probes  sets  and  heat 
map  generation  were  performed  in  GeneMaths 
XT  following  row  mean  centering  of  log2  trans¬ 
formed  MAS5.0  signal  values;  probe  set  cluster¬ 
ing  was  performed  by  the  UPGMA  method  using 
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Cosine  correlation  as  the  similarity  metric.  For 
comparative  purposes,  clustered  heat  maps 
included  probe  sets  for  spike-in  controls 
(Genisphere),  or  endogenous  small  RNAs  exhib¬ 
iting:  1)  Array  detection  p-values  <  0.05,  and  2) 
either  a)  a  log2  signal  value  standard  deviation 
<  0.025  for  all  samples  or  b)  in  the  DA  top  50%- 
tile  with  an  FC  >  1.3  in  the  opposite  direction  of 
the  selected  SDE  profile. 

Confirmatory  quantitative  Real-Time  Poly¬ 
merase  Chain  Reaction  (qPCR)  studies 

We  identified  specific  target  sequences  for  each 
of  the  18  candidate  small  RNA  biomarkers 
(Table  2).  Based  on  the  sequence  information, 
qPCR  primer  sets  specific  for  each  of  the  bio¬ 
markers  were  custom-designed  and  synthesized 
commercially  by  Applied  Biosystems  (Carlsbad, 
CA).  One  microgram  of  total  RNA  was  used  to 
prepare  complementary  DNA  (cDNA)  libraries 
using  the  High-Capacity  cDNA  Archive  Kit 
(Applied  Biosystems,  Foster  City,  California)  in  a 
total  volume  of  20  pL.  Data  were  normalized 
relative  to  those  for  the  58S  ribosomal  RNA. 
Levels  of  targeted  small  non-coding  RNA  were 
expressed  relative  to  those  in  control  groups 
using  the  2‘MCt  method  [29]. 

Results 

Using  the  Affymetrix  Fluman  gene  1.0  ST  Array 
chip  as  a  high-throughput  platform,  we  analyzed 
the  expression  profile  of  1500  small  non-coding 
RNAs  in  our  human  PBMC  specimens  from  the 
mTBI  and  non-mTBI  cases  in  our  OEF/OIF  Vet¬ 
eran  study  cohort.  We  detected  a  total  of  428 
small  RNA  species  from  our  PBMC  specimens: 
190  miRNAs,  220  snoRNAs,  8  small  cytoplas¬ 
mic  RNAs  and  10  ribosomal  RNAs.  In  an  initial 
exploratory  data  analysis,  we  conducted  a  prin¬ 
cipal  component  analysis  of  all  signals  detected 
to  assess  the  potential  value  of  the  data  sets  to 
segregate  mTBI  and  non-mTBI  cases.  The  analy¬ 
sis  revealed  that  data  sets  from  all  cases  can 
be  clustered  into  an  mTBI  or  non-mTBI  group, 
with  the  exception  of  mTBI  case  #33811,  which 
is  plotted  far  from  the  mTBI  and  non-mTBI  clus¬ 
ters  (Figure  1).  Results  from  the  principal  com¬ 
ponent  analysis  suggested  that  the  dataset  gen¬ 
erated  from  case  #33811  was  an  outlier.  Based 
on  this  and  the  fact  that  the  quality  of  RNA  ex¬ 
tracted  from  this  case  was  poor  (not  shown),  we 
excluded  case  #33811  from  subsequent  statis¬ 
tical  analyses. 
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Table  2.  Targeted  nucleic  acid  sequences  selected  for  qPCR  analysis.  Shown  are  names  of  the  18  candidate  small 
RNA  biomarkers  and  the  corresponding  targeted  nucleic  acid  sequence  used  for  the  construction  of  selective  probe 
systems  for  qPCR  studies.  Also  shown  is  a  target  sequence  for  58SrRNA,  which  we  selected  as  an  internal  control 
for  qPCR  studies. 

Small  RNA  Target  Sequence 

name 

U8  ATCGTCAGGTGGGATAATCCTTACCTGTTCCTCCTCCGGAGGGCAGATTAGAACATGAT- 

GATTGGAGATGCATGAAACGTGATTAACGTCTCTGCGTAATCAGGACTTGCAACACCCTGATTGCTCCTGTCT 
GATT 

U58B  CTGCGATGATGGCATTTCTTAGGACACCTTTGGATTAATAATGAAAACAACTACTCTCTGAGCAGC 

U27  ACT  CCAT  GAT  G  AACACAAAAT  G  ACAAG  CAT  ATGGCT  G  AACTTT  CAAGT  G  ATGT  CAT  CTT  ACTACT  G  AG  AAGT 


U83A 

HBII-289 

U55 

HBII-239 

U38B 


GCTGTTCGTTGATGAGGCTCAGAGTGAGCGCTGGGTACAGCGCCCGAATCGGACAGTGTA- 
G  AACCATT  CT  CT  ACTGCCTTCCTT  CT  G  AG  AACAGC 

ACTG AG GAATGATG ACAAG AAAAG GCCG  AATTGCAGTGTCTCCATCAGCAGTTTGCTCTC- 
CATGGGCACACG  AT  G  ACAAAAT  AT  CCT  G  AAGCG  AACCACT  AGT  CT  G  ACCT  CAGT 

GTGTATGATGACAACTCGGTAATGCTGCATACTCCCGAGTGCGCGGTGGGGAAGCCAACC- 

TTGGAGAGCTGAGC 

TGTGTGTTGGAGGATGAAAGTACGGAGTGATCCATCGGCTAAGTGTCTTGTCACAATGCT- 
G  AC  ACT  C  AAACT  G  CTG  ACAG  CACACG 

TCTCAGTGATGAAAACTTTGTCCAGTTCTGCTACTGACAGTAAGTGAAGATAAAGTGTGTCTGAGGAGA 


U56 


CCACAATGATGGCAATATTTTTCGTCAACAGCAGTTCACCTAGTGAGTGTTGAGACTCTGGGTCTGAGTGA 


U15B 

U35A 

ACA48 

U91 

ENSG0000019 

941 


CTTCAGTGATGACACGATGACGAGTCAGAAAGGTCACGTCCTGCTCTTGGTCCTTGTCAG- 

TGCCATGTTCTGTGGTGCTGTGCACGAGTTCCTTTGGCAGAAGTGTCCTATTTATTGATCGATTTAGAGGCAT 

TTGTCTGAGAAGG 

GGCAGATGATGTCCTTATCTCACGATGGTCTGCGGATGTCCCTGTGGGAATGGCGACAAT- 

GCCAATGGCTTAGCTGATGCCAGGAG 

TGTCCCTGACCTGGGTAGAGTGGCATCTGGTTGGTGATGCCCATCTCATATCAGCCAGG- 

GACAAAGCAACTCCTTGTTCATCCCAGCTTGGCTTTTGATCCGTGCCCATGCCTGGTTCATGCCTTGGACACA 

TAG 

TGGCCGATGATGACGAGACCACTGCGCAATCTGAGTTCTGGGAACCAGGTGATGGAGTAT- 
GTT  CT  G  AG  AACAG  ACT  GAG  GCCG 

TCTTAGTGACATAATTCTAATAGTTTGTTCCGACCTTCCACTGTGGACTCAATAGCAGG- 
GAG  AT  G  AAG  AG  G  ACAGT  G  ATTGCAT  G  A 


hsa-miR-671-5p  AGGAAGCCCUGGAGGGGCUGGAG 


hcmv-miR-US4  CGACAUGGACGUGCAGGGGGAU 
hsa-mi  R-1285  UCUGGGCAACAAAG  UGAGACCU 


hsa-miR-455-3p  GCAGUCCAUGGGCAUAUACAC 

58SrRNA  CGACTCTTAGCGGTGGATCACTCGGCTCGTGCGTCGATGAAGAACGCAGCTAGCTGCGA- 

GAATTAATGTGAATTGCAGGACACATTGATCATCGACACTTCGAACGCACTTGCGGCCCCGGGTTCCTCCCG 
GGGCTACGCCTGTCTGAGCG 


Using  high  throughput  small  non-coding  RNA 
datasets,  we  continued  to  search  for  candidate 
mTBI  biomarkers  that  are  differentially  regu¬ 
lated  in  PBMCs  of  mTBI  compared  to  non-mTBI 
cases.  Two  criteria  were  used  to  identify  candi¬ 
date  small  RNA  biomarkers  for  mTBI:  1)  group 
changes  (mTBI  vs.  non-mTBI  groups)  must  be 
associated  with  a  magnitude  of  >  1.5-fold,  and 
2)  group  changes  must  be  statiscally  significant 
with  p  <  0.05,  based  on  t-test  analysis  followed 
by  the  application  false  discovery  rate  correc¬ 
tions  for  multi-sampling  errors.  We  identified  18 
candidate  small  RNA  biomarkers  meeting  both 
criteria:  4  miRNAs,  13  snoRNAs  and  1  small 
scaRNA.  Interestingly,  we  observed  that  all  can¬ 


didate  small  RNA  biomarkers  are  significantly 
down-regulated  in  mTBI  versus  non-mTBI  cases. 
In  an  unsupervised  clustering  analysis  using 
RNA  expression  data  generated  from  the  high- 
throughput  gene  chip  platform  for  the  18  candi¬ 
date  small  RNA  biomarkers,  we  were  able  to 
correctly  segregate  all  17  mTBI  and  non-mTBI 
cases  analyzed  (Figure  2),  implicating  their  po¬ 
tential  value  as  surrogate  biological  indices  for 
mTBI  among  the  OEF/OIF  Veteran  population. 

Independent  qPCR  validation  of  candidate  small 
non-coding  RNA  biomarkers 

We  next  used  an  independent  quantitative  real- 
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Figure  1.  Principal  component  analysis  of  PBMC  small  RNAs  from  mTBI  and  non-mTBI  veteran  cases.  Small  RNA  ex¬ 
pression  profiles  for  individual  cases  were  assessed  by  a  high-throughput  Affymetrix  Human  gene  1.0  ST  Array  chip 
platform,  which  detected  428  small  RNAs  from  PBMC  specimens.  Signals  from  all  small  RNA  detected  for  each  of  the 
9  mTBI  and  9  non-mTBI  cases  were  summarized  into  single  points  (represented  by  balls)  plotted  on  a  3-dimensional 
plot.  Blue  and  yellow  balls  represent,  respectively,  mTBI  and  non-mTBI  cases.  The  analysis  revealed  that  all  cases  can 
be  clustered  into  an  mTBI  or  a  non-mTBI  group,  with  the  exception  of  mTBI  case  #33811  (indicated  by  a  red  arrow), 
which  is  plotted  far  away  from  the  mTBI  cluster. 


time  polymerase  chain  reaction  (qPCR)  proce¬ 
dure  to  assess  the  expression  of  individual  can¬ 
didate  small  RNA  biomarkers  in  PBMC  speci¬ 
mens  from  the  same  9  mTBI  and  9  non-mTBI 
cases  used  in  our  high-throughput  biomarker 
discovery  studies.  Primer  sets  specific  for  each 
of  the  biomarkers  were  custom-designed  and 
synthesized  commercially  by  Applied  Biosys¬ 
tems  (Carlsbad,  CA).  Using  these  primer  sets, 
we  conducted  qPCR  studies  and  assessed  the 
contents  of  individual  candidate  small  RNA  bio¬ 
markers  in  PBMCs  of  mTBI  compared  to  non- 
mTBI  cases.  Results  from  our  qPCR  studies  con¬ 
firmed  that  13  of  the  18  candidate  small  RNA 
biomarkers  are,  indeed,  differentially  regulated 
in  PBMCs  of  mTBI  compared  to  non-mTBI  vet¬ 
eran  cases  (Figure  3).  The  13  confirmed  small 
RNA  biomarkers  include  12  small  nucleolar  RNA 
(ACA48,  ENSG199411,  HBII-239,  HBII-289, 

U15B,  U27,  U35A,  U55,  U56,  U58B,  U83A,  U91) 
and  1  miRNA  (Has-miR-671-5p)  (Figure  3).  Con¬ 
sistent  with  evidence  from  our  RNA  array  plat¬ 
form,  each  of  the  13  confirmed  small  RNA  bio¬ 
markers  are  found  in  significantly  lower  levels  in 
PBMC  specimens  from  mTBI  compared  to  non- 
mTBI  veteran  cases  (Figure  3). 

Exploring  the  potential  value  of  novel  small  RNA 
biomarkers  for  segregating  mTBI  from  non-mTBI 
veteran  cases 

Based  on  results  from  our  qPCR  biomarker  con¬ 


firmation  studies,  we  next  assessed  the  value  of 
the  13  confirmed  small  RNA  biomarkers  as  a 
criterion  to  correctly  diagnose  mTBI  versus  non- 
mTBI  veteran  cases.  Using  an  unsupervised 
clustering  analysis,  we  found  the  13  confirmed 
mTBI  biomarkers  effectively  segregated  the 
cases  correctly  into  mTBI  and  non-mTBI  groups, 
with  the  exception  of  2  of  the  non-mTBI  cases 
which  were  incorrectly  identified  as  mTBI  (Figure 
4A).  We  continued  unsupervised  clustering 
analyses  to  test  the  efficacy  of  individual  or  a 
combination  of  qPCR  confirmed  biomarkers  to 
correctly  segregate  mTBI  and  non-mTBI  veteran 
cases.  Outcomes  from  these  analyses  led  to  the 
identification  of  a  3  small  nucleolar  biomarker 
panel,  comprised  of  HBII-289,  ENSG199411 
and  U35A,  which  is  capable  of  distinguishing 
mTBI  from  non-mTBI  veteran  cases  with  89% 
accuracy,  82%  selectivity  and  78%  specificity 
(Figure  4B). 

Discussion 

Evidence  has  suggested  that  appropriate  inter¬ 
ventions  can  reduce  functional  impairment  after 
mTBI  [30,31,32,33].  In  order  to  demonstrate 
the  efficacy  of  clinical  interventions,  research 
must  identify  the  biological,  clinical,  and  neuro¬ 
logical  indices  that  are  sensitive  to  the  detec¬ 
tion  of  functional  impairments  after  mTBI.  Re¬ 
sults  from  our  studies  led  to  the  identification  of 
13  novel  clinically  accessible  small  RNA  mTBI 
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Figure  2.  Unsupervised  clustering  analysis  of  18  candidate  small  RNA  TBI  biomarker  species.  The  18  differentially- 
regulated  small  RNAs  identified  from  interim  high-throughput  Array  Chip  analysis  of  8  mTBI  and  9  non-mTBI  cases 
were  subjected  to  unsupervised  hierarchical  clustering  analysis  using  the  UPGMA  algorithm  with  cosine  correlation  as 
the  similarity  metric.  Results  are  presented  as  a  heat  map  (left  panel)  demonstrating  that  the  panel  of  18  small  RNA 
biomarker  species  is  able  to  correctly  segregate  mTBI  from  non-mTBI  cases.  Names  for  each  of  the  small  RNA  bio¬ 
marker  species  are  identified  under  “Probe  Set  Name”.  Small  RNA  classes  (and  subclasses)  that  these  18  differen¬ 
tially-regulated  mTBI  biomarkers  belong  to  are  shown  under  “Probe  Set  Type”.  Vertical  dendrogram  (right  panel)  pre¬ 
sents  average  (+/-  SD)  signal  detections  from  mTBI  versus  non-mTBI  groups  for  each  of  the  18  candidate  small  RNA 
biomarkers  and  confirmed  divergent  regulations  of  the  biomarkers  in  PBMC  specimens  from  mTBI  vs.  non-mTBI 
groups.  Differential  regulations  of  the  18  candidate  biomarkers  likely  reflect  true  biological  effects  and  not  systematic 
experimental  artifact(s)  since  there  are  no  observable  group  differences  for  the  detection  of  spike-in  control  oligonu¬ 
cleotides  in  all  17  OIF/OEF  veteran  cases  analyzed  (see  heat  map  and  vertical  dendrogram).  Abbreviations:  miRNA, 
microRNA;  snoRNA,  small  nucleolar  RNA;  C/D  Box,  the  C/D  box  subclass  of  small  nucleolar  RNA;  HAc  Box,  the  HAc 
Box  subclass  of  small  nucleolar  RNA;  scaRNA,  small  Cajal  body-specific  RNA. 


biomarkers,  including  12  small  snoRNA  and  1 
miRNA.  Using  qPCR,  we  have  independently 
confirmed  that  each  of  these  biomarkers  is  sig¬ 
nificantly  down-regulated  in  PBMC  specimens 
from  mTBI  compared  to  non-mTBI  veteran 
cases.  Among  the  13  mTBI  biomarkers,  we 
demonstrated  that  a  panel  of  3  snoRNA  -  HBII- 
289,  ENSG199411  and  U35A  -  is  capable  of 
distinguishing  mTBI  from  non-mTBI  veterans 
with  89%  accuracy,  82%  selectivity  and  78% 
specificity.  Collectively,  our  evidence  suggests 
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that  additional  applications  of  the  small  RNA 
biomarkers  we  have  identified,  particularly  the 
three  biomarker  panel,  to  current  diagnostic 
criteria  may  improve  mTBI  detection  and  pro¬ 
vide  more  sensitive  outcome  measures  for  clini¬ 
cal  trials. 

PTSD  is  commonly  co-morbid  with  mTBI  in  OEF/ 
OIF  Veterans  [34,35].  We  note  that  the  majority 
of  Veterans  with  mTBI  in  our  biomarker  study 
have  co-morbid  PTSD  and  that  our  non-mTBI 
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Figure  3.  Independent  quantitative  real-time  polymerase  chain  reaction  (qPCR)  assays  confirmed  that  13  small  RNA 
TBI  biomarkers  are  differentially  regulated  in  PBMCs  of  mTBI  relative  to  non-mTBI  control  cases.  PBMC  contents  for 
each  of  thel8  candidate  small  RNA  biomarkers  identified  by  the  high-throughput  Array  Chip  platform  in  Figure  2  were 
quantitatively  assessed  using  independent  qPCR  assays.  The  same  9  mTBI  and  9  non-mTBI  cases  (Table  1)  we  used 
in  our  initial  high-throughput  biomarker  discovery  assay  were  assessed  in  this  qPCR  biomarker  confirmation  study. 
Bar  graphs  represent  mean  small  RNA  biomarker  contents  in  the  mTBI  group  relative  to  the  non-mTBI  group;  error 
bars  represent  standard  errors.  *  False  discovery  rate-corrected  P-value  <  0.05.  qPCR  confirmed  13  small  RNA  bio¬ 
marker  species  are  significantly  down-regulated  in  PBMC  of  mTBI  compared  to  non-mTBI  veteran  cases.  These  13 
confirmed  small  RNA  biomarkers  include  12  small  nucleolar  RNA  (ACA48,  ENSG199411,  HBII-239,  HBII-289,  U15B, 
U27,  U35A,  U55,  U56,  U58B,  U83A,  U91)  and  1  miRNA  (Has-miR-671-5p)  species. 


cases  are  selected  to  match  for  PTSD  diagnosis. 
Thus,  our  identified  13  small  RNA  biomarkers 
likely  represent  biological  indices  selective  for 
mTBI.  Moreover,  mTBI  cases  in  our  biomarker 
discovery  studies  were  recruited  after  an  aver¬ 
age  interval  of  3.9  years  following  their  last  de¬ 
ployment  (deployment-to-recruitment  interval: 
ranging  from  0.7  to  10.2  years,  with  a  median 
interval  of  3.4  years)  (Table  1).  Thus,  changes  in 
the  regulation  of  these  small  RNA  mTBI  bio¬ 
markers  that  we  observed  are  not  acute  mTBI 
responses,  but  likely  represent  long-term 
physiological  consequences  subsequent  to 
mTBI  experienced  during  deployment. 

The  pathological  implication  of  our  observation 
that  select  small  nucleolar  RNA  and  miRNA  are 
differentially  regulated  in  the  PBMCs  of  mTBI 
relative  to  non-mTBI  veteran  cases  is  currently 
unknown.  Small  nucleolar  RNA  and  miRNA  are 
members  of  a  family  of  non-coding  RNAs  that 
are  involved  in  many  physiological  cellular  proc¬ 
esses  and  are  also  known  to  contribute  to  mo¬ 
lecular  alterations  in  pathological  conditions 
[36].  SnoRNAs  are  short  RNA  sequences  com¬ 
prised  of  -60-220  nucleotides.  They  are  primar¬ 
ily  known  for  their  role  as  guide  molecules  for 


site-specific  methylation  and  pseudouridylation 
of  other  RNAs,  particularly  rRNA,  as  well  as  tRNA 
and  small  nuclear  RNAs.  These  chemical  altera¬ 
tions  are  required  for  proper  rRNA  processing 
and  ribosome  function  as  well  as  for  proper 
function  of  the  spliceosome  [37].  MicroRNA  are 
short  (-22  nucleotides)  RNA  sequences  that 
bind  to  complementary  sequences  on  target 
mRNA,  thereby  blocking  translation  or  promot¬ 
ing  degradation  of  target  mRNA  [37].  SnoRNA 
and  miRNA  are  expressed  in  the  brain  and  both 
classes  of  small  RNAs  have  been  implicated  in 
neuroplasticity  mechanisms  and  neurological 
disorders.  For  example,  recent  evidence  sug¬ 
gests  a  role  of  the  HBII52  small  nucleolar  RNA 
in  the  regulation  of  alternative  splicing  of  the 
serotonin  2c  receptor  [38],  and  that  patients 
with  autism  and  Prader-Willi-like  characteristics 
are  found  to  have  reduced  levels  of  HBII52  in 
the  brain  [39].  The  miRNA  miR132  is  induced 
by  neuronal  activity  and  neurotrophins  in  a 
CREB-dependent  manner  and  plays  a  role  in 
regulating  neuronal  morphology  and  cellular 
excitability  [40],  Moreover,  preclinical  evidence 
in  rodent  models  demonstrated  that  small  RNA 
expression  is  affected  in  the  brain  by  TBI.  Redell 
et  al.  [8]  reported  transient  elevated  expression 
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The  efficacy  of  the  3  Biomarker  panel  to  dissect  mTBI  from  non-TBI  veteran  cases 
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Figure  4.  The  content  of  small  RNA  biomarkers  in  clinically  accessible 
PBMCs  provides  a  sensitive  and  specific  criterion  for  dissecting  mTBI 
from  non-mTBI  veteran  cases.  We  tested  the  role  of  the  13  qPCR  con¬ 
firmed  small  RNA  mTBI  biomarkers  as  a  criterion  for  distinguishing  mTBI 
from  non-mTBI  veteran  cases.  Biomarker  contents  in  banked  PBMC 
specimens  from  the  same  9  mTBI  and  9  non-mTBI  cases  (Table  1)  were 
used  in  our  biomarker  discovery  studies  and  quantified  by  qPCR.  (A,B) 
The  efficacy  of  using  biomarker  contents  from  clinically  accessible 
PBMCs  as  a  criterion  to  correctly  segregate  mTBI  and  non-mTBI  cases 
was  tested  by  unsupervised  clustering  analysis  using  the  UPGMA  algo¬ 
rithm  with  cosine  correlation  as  the  similarity  metric.  Results  are  pre¬ 
sented  as  heat  maps  demonstrating  the  efficacy  of  using  all  13  small 
RNA  biomarkers  (A)  or  using  a  panel  of  three  small  nucleolar  RNA  bio¬ 
markers  (B)  to  correctly  segregate  mTBI  from  non-mTBI  cases.  (B)  A  three 
small  nucleolar  RNA  biomarker  panel  (HBII-289,  ENSG199411  and 
U35A)  is  capable  of  distinguishing  mTBI  from  non-mTBI  cases  with  89% 
accuracy,  82%  selectivity  and  78%  specificity  (accuracy  is  the  percentage 
of  all  mTBI  and  non-mTBI  cases  that  are  correctly  identified;  sensitivity  is 
the  probability  that  a  case  identified  as  mTBI  actually  is  a  mTBI  case; 
specificity  is  the  probability  that  a  case  identified  as  a  non-mTBI  case  is 
actually  a  non-mTBI  case).  Abbreviations:  snoRNA,  small  nucleolar  RNA; 
miRNA,  microRNA. 


miRNA  species  in  the  rat  cere¬ 
bral  cortex  up  to  72  hrs  following 
fluid  percussion  injury  to  the 
brain.  Redell  et  al  [42]  also  iden¬ 
tified  potential  up-regulation  of 
35  and  down-regulation  of  50 
miRNA  species  in  the  hippocam¬ 
pus  of  rats  within  72  hrs  follow¬ 
ing  an  impact  TBI  to  the  brain; 
altered  regulations  for  a  smaller 
subset  of  8  (4  up-regulated  and 
4  down-regulated)  miRNA  spe¬ 
cies  in  the  hippocampus  were 
subsequently  confirmed  by  inde¬ 
pendent  qPCR. 

It  is  possible  that  the  altered 
regulation  of  select  small  nucleo¬ 
lar  RNA  and  miRNA  that  we  ob¬ 
served  in  PBMCs  of  veteran 
mTBI  cases  might  have  implica¬ 
tions  in  the  central  nervous  sys¬ 
tem.  Genes  relevant  to  neural 
circuits,  synapses  and  neural 
plasticity  processes  are  also  ex¬ 
pressed  in  circulating  blood 
cells,  such  as  PBMCs.  Thus,  the 
fact  that  we  observed  significant 
down-regulation  of  select  small 
RNA  biomarkers  in  PBMC  speci¬ 
mens  from  our  veteran  mTBI 
cases  long  after  their  deploy¬ 
ment  might  reflect  long-term 
molecular  alterations  in  the  cen¬ 
tral  nervous  system  contributing 
to  the  onset  and  progression  of 
clinical  TBI  phenotypes.  Future 
studies  exploring  the  physiologi¬ 
cal  relevance  of  mTBI  bio¬ 
markers  will  provide  a  better 
understanding  of  the  biological 
mechanisms  underlying  mTBI 
and  insights  into  novel  therapeu¬ 
tic  targets  for  mTBI. 
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